Has the Universe always expanded ? 
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We consider a cosmological setting for which the currently expanding era is preceded by a con- 
tracting phase, that is, we assume the Universe experienced at least one bounce. We show that 
scalar hydrodynamic perturbations lead to a singular behavior of the Bardeen potential and/or its 
derivatives (i.e. the curvature) for whatever Universe model for which the last bounce epoch can be 
smoothly and causally joined to the radiation dominated era. Such a Universe would be filled with 
non-linear perturbations long before nucleosynthesis, and would thus be incompatible with obser- 
vations. We therefore conclude that no observable bounce could possibly have taken place in the 
early universe if Einstein gravity together with hydrodynamical fluids is to describe its evolution, 
and hence, under these conditions, that the Universe has always expanded. 

98.80.Hw, 98.80.Cq 



I. INTRODUCTION 

Since the development of singularity theorems 0|, it 
has been a paradigm that the Standard Big Bang (SBB) 
cosmology [Q| requires the Universe to have emerged from 
an initial singularity. However, as all theorems, those 
concerning singularities were based on a set of hypoth- 
esis. In the 1970's and after, many non-singular_cosmo- 
logical models with bounces were constructed [|| Jj where 
indeed one or more of the hypothesis of the theorems 
were violated, such as, for instance, energy conditions or 
the validity of Einstein gravity. 

More recently, superstring-inspired models suggested 
that this singularity might have been avoided if gravity 
was modified somehow in the near Planckian regime for 
which a maximal curvature was supposed to be reach- 
able. Note that although this model seems to have al- 
ways expanded in the Jordan's frame, a bounce behav- 
ior is also present in the Einstein frame. Such a pre 
Big Bang (PBB) scenario ^ is, however, plagued with 
various problem, as, e.g. the graceful exit and fails 
to reproduce the primordial perturbation spectrum as 
observed in the recent Cosmic Microwave Background 
(CMB) data §. The non-singular hypothesis was revig- 
orated even more recently with the Ekpyrotic Universe 
proposal 1^. 

Quantum gravity also provides new ways out of the sin- 
gularity theorems. For instance, models exist in which an 
expanding Universe region might originate from a Black 
Hole collapse jo). Quantum cosmology, by the Wheeler- 
de Witt equation, also yields bouncing behaviors for the 
scale factor [ p^lpd] ]. However, in this case, whether the 
Universe actually avoids the singularity may depend on 
the interpretation that is given to its wave function. In 



*Already in the 1930's, 
was proposed by Tolman 
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particular, in the many- world view p^ , it is the expec- 
tation value of the scale factor that bounces, so that, 
in practice, nothing can be said as to the value it takes 
in the realization corresponding to our Universe. How- 
ever, in the ontological interpretation of the wave 
function of the Universe, the Bohmian trajectory, being 
meaningful, implies an actual physical bounce. (Note the 
interesting fact that this interpretation also offers a nat- 
ural solution to the graceful exit problem in the PBB 
scenario.) 

In this paper, we restrict our attention to pure gen- 
eral relativity (GR) and examine a Friedman-Robertson- 
Walker (FRW) setting whose stress-energy source comes 
from a collection of perfect fluids (having vanishing 
anisotropic stresses). We then assume that a bounce took 
place at some instant in time and concentrate on observa- 
tionally meaningful models, namely those for which the 
bounce epoch can be smoothly and causally connected 
to a radiation dominated phase where Big Bang Nucle- 
osynthesis (BBN), assumed as a physical requirement in 
whatever model, can safely take place. Such Universes 
will be called "realistic" hereafter. What precedes the 
bounce itself is also left out of our analysis which is thus 
not restricted to non-singular models. In short, we just 
assume that, for some reason, our Universe has not al- 
ways expanded, and we question the observable conse- 
quences of this possibility. 

This work is to be contrasted with Ref. ||l^] in which 
the mathematical conditions on the total fluid for the 
bounce to occur are derived (see the Appendix for a 
physical restriction of these conditions) : our analysis 
is in some sense a follow-up of this work as we examine 
stability. 

In the framework defined above, we investigate the be- 
havior of scalar adiabatic perturbations near the bounce 
through the Bardeen gauge-invariant potential. We show 
that these perturbations can grow without limit in two 
situations. One, when the Universe passes through the 
bounce itself. The other, if the existence of the bounce 
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imposes a violation of the Null Energy Condition (NEC) 
(and hence of all the other energy conditions ||l5| ) , when 
the Universe makes the transition from the regiofl where 
the NEC is violated to the region where it is notJll. In all 
cases of physical interest, at least one of these instabilities 
occurs. These facts, in turn, imply that there must exist 
a point in time where, apart for infinitely fine-tuned ini- 
tial conditions (and yet for all wavenumbers at the same 
time), the perturbations grow unbounded. This means 
that if such a bounce had occurred, then very large inho- 
mogeneities, at all scales, would have been formed very 
early on, that would presumably have turn non linear 
very rapidly, i.e. long before nucleosynthesis. Such a 
conclusion is in flagrant contradiction with BBN, so that 
we are led to postulate a no-bounce conjecture ||l8| . 

In fact, our conclusion is more general than that. Even 
without bounce, the point (in time) where the NEC is re- 
covered (assuming it has been violated in the past, which 
is the case for the majority of bouncing Universes) yields 
divergences in the primordial perturbations. As a result, 
the PBB models, which may have no bounce, acquires a 
new constraint as follows : the graceful exit problem, in 
order to be solved, also requires, in many instances p9[ , 
a violation of the NEC for a finite amount of time (see 
however Ref. Q ) . Models should then take into account 
the fact that the recovery of the NEC should not occur 
at a time where our hypothesis (GR plus hydrodynam- 
ical perturbations) are already valid to avoid the excess 
perturbation problem. 

Finally, there is the special case, which is actually not 
a bounce, for which the scale factor, although never de- 
creasing in time, happens to stop at an instant (a(?7o) = 
ao + bijl^^^ + ... for some 770, with i > 1). Such a "pseudo- 
bounce" situation is included in our analysis, and it 
also generates a cosmological catastrophe. One must 
therefore conclude that the expansion must be strictly 
monotonous in order to make sense. 

In what follows, after a short review (Sec. 0) on the es- 
sentials of cosmological perturbation theory in the frame- 
work with which we are concerned, heavily based on 
Ref. M, we detail the structure of a bouncing epoch 
(Sec. [II) by making an expansion of the scale factor in 
powers of the conformal time. We then go on to eval- 
uate the evolution of the primordial perturbations (in a 
gauge invariant manner in order to avoid spurious modes) 
around the bounce in Sec. establishing our results. 
Sees. and Wll are devoted to discussions and conclu- 
sions, while a final appendix shows that connecting the 
bouncing phase with the observable Universe yields con- 
straints on the parameters appearing in the expansion 



^Fluids which violate the NEC are called exotic matter. 
They have strange properties but there is no theorem assert- 
ing their impossibility on macroscopic scales jig] (one way or 
the other) For examples of energy-momentum tensors where 
the NEC is violated, see Ref. U7\ 



and exhibits clearly the violation of the NEC in many 
cases. 



II. CLASSICAL COSMOLOGICAL 
PERTURBATION THEORY AND THE NULL 
ENERGY CONDITION 

The purpose of this section is to recall various defi- 
nitions and evolution equations for scalar hydrodynamic 
perturbation, for the sake of self-consistency. 

Cosmological perturbations are expected to be pro- 
duced quantum mechanically during the early stage of 
the evolution of the Universe. Expansion drives the quan- 
tum mechanical fluctuations of any field present in one's 
favorite model to classical values, which then evolve ac- 
cording to GR. Although we will concentrate on the sub- 
sequent evolution of primordial perturbations as soon as 
they have been produced, we shall also summarize the ac- 
tual process of building them as the quantum modes of 
the field. This quantum analysis yield different insights 
into the nature of the problem. 

The background metric we shall be using is 



ds^ = 0^(77) (dry^ - 7ydxMx^ ) , 



(1) 



with rj the conformal time, related to the cosmic time t 
through the scale factor a, adi] = dt, and 



7y 



(2) 



the 3-space metric. Here, K, = 0, ±1 is the sign of the 
curvature. The general form of metric perturbations on 
this background reads 



ds^ = a2(77) {(1 + 2(j))d7j^ - 2B,idr]dx' 
- [(1 - 2i/.)7y + 2E,,j] dx'dx^} 



(3) 



From this form, it is convenient to define the gauge- 
invariant Bardeen potentials pl| 



^[iB-E')a]', 



4' = V---(S-S'), 
a 



(4) 



(5) 



(where a prime means a derivative with respect to con- 
formal time), in terms of which it is possible to expand 
the perturbed Einstein equation in a gauge independent 
way. The latter are further simplified in the particular 
case of vanishing anisotropic stress for which 5* = <I>. In 
what follows, we shall restrict our attention to this spe- 
cial case. Note that for perturbation theory to be valid, 
and space not to become very inhomogeneous, the re- 
quirement |<I>| ^ 1 should hold at all times and at all 
wavelengths of perturbation. The stress energy tensor, 
source of Einstein field equations, will take the form 
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(6) 



for energy density e, pressure p and 4- velocity u. Follow- 
ing Ref. Q|, one then defines the Null Energy Condition 
through the Null Convergence Condition which requires 
the Ricci tensor to satisfy R^un^^n" > for all null vector 
. If this condition is fulfilled, then Einstein equations 
imply that, with the form (|6|) of the stress-energy tensor, 
the relation 



NEC 



e+p > 0, 



(7) 



holds. This relation, whose violation implies the violation 
of all other energy conditions, is indeed assumed to hold 
in almost every situation of physical relevance. As we 
shall see, a bouncing Universe requires that it be violated 
in many circumstances. 

Perturbations of the stress energy tensor read 



_ eo + Po 



Sui 



ST' = -(5p(5". (8) 



Here, and in the rest of this paper (unless otherwise 
stated), an index on a quantity is meant to be the back- 
ground, unperturbed value of this quantity. These per- 
turbations are made gauge-invariant through the trans- 
formation 6e Jets') = Se + e'„{B- E'), 5p Jpfe') = 

Sp + p'o{B - E'), and 5u^ 5uf'^ = 5u, -f a{B - E').,. 
Einstein equations then read, in terms of the Bardeen 
potential 



(a$)', = 47rGa^(eo + po)Su^ 



(gi) 



(9) 



(10) 



$" + 3n<i>' + {2-H' + H^- /C)$ - AirGa^Sp^s'l (11) 
where the background equations (Ti. = a' /a) 
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+ IC — —TrGa^eo, 



-P-n'-JC = -7rGa'po, 



(12) 



(13) 



have been taken into account. In Eq. (|l^), we have used 
the variable f3 defined as 



(14) 



which we will use intensively later on. It can be remarked, 
at this stage, that Eqs. (03) and (HSh imply 



eo + Po = 



/3 



47rGa2 ' 



(15) 



so that a violation of the Null Energy Condition Eq. (|^) 
will be equivalent to having f3 <0. 



In general, pressure perturbations can be separated 
into adiabatic and entropic perturbations through 

with 6S the entropy perturbation. In many situations, T 
is interpreted as the squared sound velocity of the fluid 
and denoted . However, for the fluids which appear in 
the examples below, this variable cannot be interpreted 
as such. Already in plasma physics there are situations 
where T^/^ is not the group velocity of the fluid (see 
Rcf. page 406, for a discussion on this point, and 
references therein). That is why we shall avoid calling 
-s/T the sound velocity in what follows. 

Using Eq. (|l|), and plugging Eq. (0) into (|), one gets 

+ 37i:(l + T)$' - TV^* 
+ [2H' + (1 + 3T){H^ - /C)] $ = 47rGaV(55, (17) 

which becomes an homogeneous equation in the case of 
adiabatic perturbations ^.S* = with which we shall be 
conc erne d from now on (see however the discussion in 
Sec.[VB|). 

We shall need, in what follows, two more pieces of in- 
formation. It can be seen, by combining Eqs. ([l^), ( [T^ ) 
and (HQ) that the following useful identity for T 



/3' 



(18) 



holds. Moreover, it is clear that perturbation theory can 
only make sense if the perturbed Einstein tensor remains 
a small quantity. It turns out that Einstein equations do 
not involve anything but the Bardeen potentials. This 
stems from the fact that the gauge-invariant part of the 
Einstein tensor SG^f^^'^ [the left-hand side of Eqs. (||) to 
(prj)] only depends on $ and ^' |Q. Satisfying the gen- 
eral requirements for perturbative analysis thus demands 
that not only the Bardeen potential, as discussed above, 
but also its first and second derivative in time, in order 
to remain a perturbation, should be regular and well- 
behaved for all times. As we shall see later, this is not 
always a trivial requirement. 

Hydrodynamical perturbations can be shown to origi- 
nate from a single gauge-invariant variable v{ri, x), whose 
dynamics is derivable from the simple action 



S" = — y ^y^d'^x ^v'^ — TY'' v^iVj H 1 



(19) 



where 7 is the determinant of the 3-metric 7^ (with in- 
verse 7*-'), and the background function z{ri) is 



(20) 



with f3 defined in Eq. (14). The variable v is formed with 
the fluid velocity perturbation potential and the Bardeen 
potential, its explicit form being here however irrelevant. 
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Eq. ( |l9| ) is in the correct form for a straightforward 
canonical quantization of the field v. There exists a con- 
nection between v and the Bardeen potential which we 
will use below and that we therefore reproduce here. We 
first expand $ on a basis of eigenfunctions of the Laplace- 
Bcltrami operator (here, fc is a generic eigenvalue of this 
operator, which, depending on the curvature /C, can be 
continuous or discrete). 



I dk {M^)nl{fl)au + h.c.) , (21) 



where (. = a/SttG/S is the Planck length, and V'/c(x) sat- 
isfies 



(V2+fc2)^fe(x)=0. 

Implementing the usual commutation relations 



^kk' 



(22) 



(23) 



(the 5 here being a Dirac distribution or a Kronecker 
symbol depending on whether the space is flat or open, 
or closed, respectively) yields the time evolution equa- 
tion for the wavelength-dependent modes of the Bardeen 
potential 



e 



ukiri) = 0, e 



(24) 



An expansion similar to Eq. ( pi|) can be done for the field 
V. It turns out that its modes Vk{vi) satisfy 



4'(^)+(Tfc2--)«fe(^)-0, 



(25) 



and the^ are related to the Bardeen operator through the 
relatior 



a\zJ 



(26) 



A last interesting point concerning this field is that the 
eigenmodes are expected to be normalized as 



K{v)v*{v) - Vk{v)vk{v) = 2i, 



(27) 



for all times rj. 

The process of generating primordial fluctuations then 
goes on by requiring the Universe to be empty at some 
early time, i.e. to be in the vacuum state |0) for which 
afc|0) = for all wavenumber k. The evolution of v 
through Eq. ( p5| ) then necessarily changes the state and 
dynamically produces particles that are subsequently in- 
terpreted as primordial perturbations. 

The background equations presented in this section 
will now be used for the forthcoming analysis of the 
bouncing background, to which we now turn. 



"''When /3 or T change sign , it is necessary to calculate $ 
directly by means of Eq. (|l7[) . 



III. GENERAL BOUNCING BEHAVIOR 

In this section we will enumerate some properties of 
a background spacetime with a bounce, and calculate 
some relevant quantities around the bounce which will 
be important for the calculations of the perturbations. 
We shall be mostly concerned in what follows by actual 
bounces, but the pseudo-bounce situation is also included 
in our description. 

Our sole assumption on the background is that, near 
the last bounce before the standard model evolution took 
place, a(ry) can be written as follows 



a = flo + 677^" + d77^"+^ + eyy 



„2n+2 



(28) 



where the integeil n satisfies n > 1. It means that a(r]) 
must be at least 0^"+^ near the bounce. In order that 
Eq. ( p8| ) indeed represents a bounce, the otherwise arbi- 
trary parameter & > 0. Our results are independent on 
what happens before this last bounce (other bounces, a 
PBB evolution, or even a singularity). In the case where 
eo — eo(a)j which is possible if and only if a{r]) is even, 
the presence of a bounce implies that the model is free of 
singularities. This is the case of models of Refs. 

For this model to be realistic, we require oq to be less 
(or even much less) than a„ ^ 10^^ cm, where a„ is the 
value of the scale factor in the beginning of nucleosyn- 
thesis. 

The function Ti.{ri) coming from Eq. (|2^ ) reads 



7i = — 2 [2nbaori 



2n-l 



(2n + l)aod7] 



2n 



+ {2n + 2)eao772"+i - 2n&2,y4«-i 
while P(r]) is 
1 



(29) 



;[ICaQ — ao2n(2n -|- l)dri 



2n-l 



ao2n{2n — 1)67] 



2n-2 



-aoi2n + 2){2n + l)erf'' + 2n{<6n - l)h^r]^''-% (30) 

From these functions we can calculate the coefficient T 
given in Eq. (p!8|), yielding : 

^_ ao r-ao277(27i- l)(27i-2)6 
Qnh \ rf 
+2[lCb - {2n + 2)(2n + l)e]ao7i 
2ao7i(4772 - l)d 



V 



+166V (277 -1)77^"-^ + 



2n-2 



(/Cag - aQ2n{2n - 1)677^""^ + ...) 



(31) 



We will now evaluate the behavior of the background 
functions Ti, /3, T and z near the bounce for the possible 
values the free parameters can assume. 



^Or half- integer in the case of a pseudo-bounce. 
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1. n > 1 and /C 7^ 0. 

In this case, we find the following behaviors 



T = 



+ {JCb - 12e)ao 



T = 



(27i-l)(2n-2) 



2nb 

n = — 7 

ao 

P = IC, 

1 



(32) 

(33) 
(34) 
(35) 



so that it turns out to be possible to have a bounce 
without violation of the NEC only if /C = 1. Note 
that this well as the following one, can 

accommodate a pseudo-bounce (where 2n is an odd 
number strictly greater than 1). 

2. n > 1 and /C = 0. 



For the special case of a flat background, the vari- 
ous quantities needed to describe perturbations are 
modified as 



^_ ao(2n-2) 
ao 

(3 = - — 2n(2n~ l)r;2 
ao 

z oc const. 



(36) 

(37) 

(38) 
(39) 



Here there is always violation of the NEC as /? < 
at least in some open neighborhood of the bounce. 

3. n=l,d^Q, V/C. 



This is the case where the second derivative of a{rj) 
is non vanishing and 0(77) is not even. The relevant 
quantities are 



T = 



aad 



b{2b - ICao)r]' 

n = —T], 

ao 

ao 

1 

Z GC . 



(40) 
(41) 
(42) 
(43) 



In this case, we can also have a bounce without vi- 
olation of the NEC if /C = 1 and 26/ao < 1. How- 
ever, this condition is not physically plausible, as it 
is argued in the Appendix. 

4. n = 1, V/C, and d = 0. 



Finally, for the case that can represent an even 0(77), 
we get 



36(26 -/Cao) 

n = — ?7, 

ao 

ao 

1 

z oc — . 
V 



(44) 
(45) 
(46) 
(47) 



Here again, we can have a bounce without violation 
of the NEC if /C = 1 and 26/ao < 1. 

This case corresponds to models of Refs. P,p^. 

Note that there are exceptional cases when JC = 1 and 
26 ao in cases 3 and 4, and when 86^ + (/C6— 12e)ao = 
in case 4. For case 3 with /C = 1 and 26 = ao, it can 
be shown that unbounded growth occurs in the bounce, 
but we will not exhibit these calculations below because, 
as we have already mentioned, /C = 1 and 26 = ao are 
not together plausible in realistic models. For case 4, 
we will show in the next section that the crucial instant 
of instability is not in the bounce itself, but when the 
Universe makes the transition from the region where the 
NEC is violated to the region where it is not. The oc- 
currence of this transition is independent on the value of 
862 ^_ (/C6- 12e)ao. 



IV. NON-LINEAR COLLAPSE OF NON 
SINGULAR MODELS 

Having set the general conditions under which a 
bounce is expected to occur in the early Universe, we 
now apply the cosmological perturbation theory recalled 
in Sec. l| to such a background. In particular, we shall 
make extensive use of Eqs. (^7|) and (|24-26) to derive 
the leading terms in the behavior of the gauge-invariant 
Bardeen perturbation potential. Let us turn back to the 
time the bounce is supposed to have taken place, and 
consider the various situations for the free parameters 
already discussed above. 



A. The bounce 

There are four distinct cases near the bounce which 
deserve to be examined. 

1. n > 1 and /C ^ 0. 

The first case is > 1 with non vanishing curvature 
K = ±1. Eqs. (H) to (H) transform Eq. (||) into 



„ (2n- l)(2n- 2) / fc^ 



1 Ufc = 0, (48) 
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whose solution is f = v+rf+ + v-rf- {v+{k) and 



. (fc) being two integration constants), with 2pj^ 



1 ± 



4(2n-l)(2n-2) 



- 1 



For nega- 



tive curvature space, /C = — 1, one has p- < 0, for 
whatever value of fc, while for /C = 1, p_ < pro- 
vided < 3. Note that Vfc, (p_ — 1) has a negative 
real part in both spaces IC — ±1. 

Inserting this form into Eq. (|6|) yields, for the 
Bardcen potential, $fc oc rjP- , i.e. also a divergent 
quantity for p_ < 0. It can be seen from the previ- 
ous considerations that, irrespective of the value of 
the wavenumber fc, there will be a divergent quan- 
tity either in the metric perturbation directly for 
/C = —1, or in its derivative (oc r/P-^^) for /C = 1. 

2. n > 1 and /C = 0. 



The second case, similar to the first {n > 1) but 
for flat spatial sections (/C = 0), has z = const. 
In order to calculate , the "pump" term, we need 
to obtain z in the next to leading order of approx- 
imation. However, making use of the expansion 
formula for Vk 



Vk ^ z 



rio ^ 



" d^ 

Vo z 



one can verify that Vk contains a term proportional 
fc^/yy^""^, which diverges in the bounce. 

In order to calculate what happens to the Bardeen 
potential ^k, it is easier to investigate directly the 
Bardeen equation (^7|) when these values of the pa- 
rameters are assumed, and one finds 



,, 2 - 2n , 



00(71 — l)fc^ 



whose solution, setting 



(49) 



(50) 



is given by (see Eq. (8.491/3) in Ref. ||) 

= j^(2"-i)/2 |$_^/2^(a) + $_ii:_2^(a)| , 

(51) 

with two arbitrary constants of integration (de- 
pending on k) and In this paper the sym- 
bols J, N and /, K stands for pairs of linearly in- 
dependent solutions of the Bessel equation, Bessel 
and modified Bessel functions of first and second 
kind. As, by assumption, n > 1, we find that the 



leading order behavior of the Bardeen potential is 
[Ref. ||], Eqs. (8.451/3-4)] 



^(3n-2)/2g|a|^ 



(52) 



namely, this perturbation exhibits a singular be- 
havior near the bounce (|a| — > 00). 

3. n^l, d^O,yiC. 

This is the most general case, where a(ry) is not even 
and its second derivative is not null. Eqs. (|4^-|47|) 
when inserted into Eq. (Eq), yield 



/ A(fc) _3_ 



vk = 0. 



(53) 



The explicit form of A is not important. The solu- 
tion of this equation reads [Ref. [|^, Eqs. (8. 404/1- 
2)] 

Vk^ ^[v+J2{2^)+V-N2{2^)Y (54) 

where , w_ are arbitrary constants depending on 
k. The leading order behavior of the divergent 
term is Vk (x. 1/y^. The Bardeen potential can be 
obtained from Eq. (|6|), or directly from Eq. ( |T7| ) 
(with SS = 0), yielding, near 77 = 0, 



k^^[v+Jii2^)+v-Ni{2^)} , (55) 



which is convergent at 77 = 0. However, the first 
derivative of $fc has a divergent term, 



(56) 



cx ^Noi2^), 



which behaves as Inry near 77 = 0, and the second 
derivative is 



1 

77 



(57) 



yielding unbounded growth of the perturbed Ein- 
stein tensor (curvature) near the bounce. 

4. 71 = 1, V/C, and d = 0. 



This is the case which contains models from 
Refs. ^,|l^, namely, that for which the first term 
in the scale factor expansion is quadratic {n = 1), 
for whatever value of the spatial curvature /C, and 
which can accommodate the case for which 0(77) is 



Eqs. ( 



inserted into Eq. (^5|), yield 
„ 2 



77^ 



TVk = 0, 



(58) 



yielding the divergent mode Vk oc I/7;. For the 
Bardeen potential, we can use Eq. (p4|), or directly 
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Eq. ( p7| ) (with 5S = 0), to obtain that, near = 0, 
it behaves as 



oc cos(Q!(fc, IC)r]) , 



(59) 



which is regular up to second order derivatives in- 
dependently on the form of the function a{k,lC). 

This is a case where there is no divergence in the 
bounce itself. However, as shown in the Appendix, 
near the bounce of a realistic Universe of this type, 
/? < 0, even when /C = 1. As the Universe must 
have reached a region in which the NEC must be 
satisfied at later times, there must exist an instant 
ryo at which the function (3 vanishes. As it turns 
out, this is a critical point also from the point of 
view of perturbation theory. 



B. NEC transition 

In what follows, we now concentrate on the point where 
/? = 0, so that we shift the origin of time : for now on, 
ryo = for /Jq, and we denote by an index quantities 
evaluated at this point. 

We now assume that the scale factor around 77 = is, 
again, derivable at least up to third order, so that the 
following expansion 



a{rf) = ao 



(60) 



holds. In this relation, ay, = a"'(0)/ao, and use has been 
made of 



27^0 + /C, 



(61) 



which is a simple rewriting of /?o = 0. 
Using the expansion (^), we find that 

H = Ho+ f]{Hl + JC) 



and 



(3 = {ml + 5/CHo - 03)?? + 0{Tf). 



(62) 



(63) 



As shown in the Appendix, /3q 7^ for realistic models 
and hence QTLq + 5ICHo — 03 7^ 0. As a result. 



T = 



1 



which transforms Eq. (|l^) near = into 



drf' rj drj 



Ho 'STCo 



(64) 



0, (65) 



whose solution reads ( ||2|l, Eq.(8.494/5)) 
<I>fc = $,77/2 



/ K- 



(66) 



where $j and are arbitrary constants depending on 
fc, while I2 and K2 are ordinary Bessel functions. For 
realistic models, Ti^ ^1. To lowest order, the Bardeen 
potential is then found as 



$fc ~ const. + 0(77), 



(67) 



which would seem to imply a safe perturbative expansion. 
However, its second derivative with respect to conformal 
time exhibits, irrespective of the value of fc, the mildly 
diverging behavior 



$','^ln77 



(68) 



so that, by virtue of Eqs. (^) - (pT[), we are forced to 
conclude that even this perturbative expansion generates 
unbounded growth. 



V. DISCUSSION 

The various cases that have been investigated here rep- 
resent an exhaustive list of all the possibilities that sat- 
isfy the requirements oi i) having a non-strictly mono- 
tonic behavior for the scale factor at some instant of 
time, ii) that this instant of time be in our past light 
cone, and that Hi) therefore, the Universe in which the 
bounce took place is ours, in the sense that nucleosyn- 
thesis as we know it occurred. In all these diver- 
gent behavior for the cosmological hydrodynamical per- 
turbations have been found, which are responsible for an 
observationally excluded non-linear growth of primordial 
inhomogeneities. We would like to discuss three points 
regarding these results. 



A. T ^ 00. 

The first important point we should like to emphasize 
concerns the divergent behavior for T, which could be 
expected to be responsible, right from the outset, for the 
divergence of the perturbations. That this is not the 
case can be seen on the following example, based on the 
last case of the previous section. If j3' — 0, i.e., if /? cx 
^2n+i ^[^]^ n > 1, when (3 = 0, then, to leading order, T 
behaves as T ~ — 3^^^ , and Eq. ( |65| ) is transformed in 
such a way that the Bardeen potential now reads 
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+2 



, (69) 



which has convergent derivatives up to order 2n + 2. This 
is a case where, even in the point where T diverges, the 
curvature and all other relevant quantities converge. It 
happens because 5e goes to zero at this point. Hence, 
stable models with a divergent T at one point are indeed 
possible, as this simple example testifies. It is the combi- 
nation of all the conditions i), ii) and Hi), that gives rise 
to unbounded growth. 



B. Gauge invariance 

In each of the cases investigated above, we have found 
a growing mode for the Bardeen potential, which leads 
to divergences either in the metric perturbation or in the 
curvature. In view of recent discussions of similar modes 
found in the PBB case for instance |^^, a natural ques- 
tion to ask then is whether they are actually physical. 
Here, we should like to argue that indeed, and contrary 
to their counterpart in PBB, they cannot be "gauged- 
down" . In this sense, they represent real physical inho- 
mogeneities, doomed to grow non-linear. 

It is well-known that for wavelength larger than H~^, a 
residual gauge artifact might exist that can be tamed for 
A; <C 7i (for instance by going to an off-diagonal gauge in 
which higher order derivatives are involved [Q). How- 
ever, the divergences we have obtained in the last section 
are present for all wavelengths. Furthermore, contrary to 
the PBB case, our divergences occurs on regular points of 
spacetime, where Ti. is regular (even zero in the bounce, 
implying that all wavelengths are smaller then the Hub- 
ble radius at this time), and it is not possible to blame a 
bad behavior of Ti for their existence. 

Let us see that explicitly. Assume one makes the cal- 
culations in the longitudinal gauge for which B = E — 
(note that this is not the case of the present paper). In 
this gauge, Eq. (||) simplifies to 



a^i-q) [(1 + 2(t))dri^ - (1 - 2V')7ijda;Ma;^] , (70) 



where both and (p take the gauge-invariant value $. 
We now assume that there exists a gauge transformation 
ai" = a;" + (x, 77) that could cancel the diver- 
gence in the potential leading to (in Fourier space, 
^^(77) being the /C— dependent Fourier transform of the 
infinitesimal vector field ^") = - rL£,l[ri) ~ C2'(?7), 

^fe = $fc + mt Bk = il- i'k and Ek = -?fe. The sec- 
ond of these relations reveals that since is a divergent 
quantity at the physical point of space-time under con- 
sideration, the relevant transformation must be singular 



in order to produce a bounded value for . This is to be 
contrasted with the PBB case for which, as r\pgg 0, 



77 so that the gauge vector ^ can indeed as- 



sume infinitesimal values while canceling the divergence. 
In our case, 7i — > Tio, a constant which, at the bounce, 
actually vanishes, thereby rendering the transformation 
even more singular. However we look at these perturba- 
tions, they do grow infinitely, leading to singularities, or 
at least to large (non-linear) inhomogeneities. 

This can also be seen by using the gauge invariant cur- 
vature perturbation variable p4] 



C = i> 



Hp 



(71) 



which also exhibits divergences in the cases of interest 
here as can readily be checked. However, one should use 
this variable cautiously as the corresponding divergences 
may not be physical : in deriving this form for the three- 
curvature modes, it is explicitly assumed that the NEC is 
satisfied, so that for /3 — > 0, a fake divergence can occur. 
This is indeed what happens in t he situation for which 

in this case, 



/3 oc Ty^i-i-i mentioned in Sec. VA 



c 



even though nothing particular happens in the theory at 
this point. 



C. Adiabatic and entropy perturbations 

The analysis that is presented here is restricted to 
purely adiabatic perturbations, which might be seen as 
too restrictive an hypothesis in view of the fact that at 
least two barotropic fluids are needed in a hydrodynami- 
cal description to yield a bounce. Bardeen equation ( p^ ) 
for the potential $ can be read as an inhomogeneous 
equation with a source term proportional to the entropy 
perturbation, if the latter is independent of the former. 
In this situation, its solution is given by a particular solu- 
tion of the full inhomogeneous equation, whatever it may 
be, plu s the general solution of the homogeneous one, i.e. 
Eq (n^ with (55 = 0, which is the equation we have been 
discussing throughout. As it was found that the adia- 
batic perturbations do exhibit a singular behavior near 
the bounce or at the time where the NEC is recovered, 
all the perturbations will exhibit such a behavior. 

What if the entropy perturbation depends on the 
Bardeen potential ? In that case, it is useful to recall 
that an arbitrary perturbation can be decomposed over 
its possible modes (in the linear regime with which we 
are concerned), namely one adiabatic, and — 1 isocur- 
vaturc modes for the case of A^ constituents. All these 
modes may be well behaved, but one sees that with the 
decomposition 25 



N-l 



"J* adiabatic + *^ isocurvaturc ' 



(72) 
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it is sufficient that the adiabatic mode we have been 
studying diverges to yield a cosmological catastrophe, un- 
less the isocurvature modes somehow compensate exactly 
this divergence. Here, again, we return to recent obser- 
vations of the CMB 1^ according to which the dominant 
part of the primordial perturbations was in the adiabatic 
mode. Therefore, once again, the requirement that the 
bouncing Universe be ours implies a stringent constraint : 
in this case, it is that the adiabatic mode under scrutiny 
here is not negligible, and therefore cannot be compen- 
sated by an isocurvature mode. 

VI. CONCLUSIONS 

The simplest way to envisage a cosmological model 
without an initial singularity is to think that the scale 
factor never shrunk to zero but reached a minimal value, 
preceded by a contraction phase. In fact, much of non- 
singular cosmological models proposed in the literature 
present this bouncing behavior. 

We have shown in the present paper that any bounc- 
ing solution whose matter content is described by hy- 
drodynamical fluids is highly unstable, at least if they 
are to be smoothly joined with our observed Universe. 
In particular, applying to these solutions the now stan- 
dard cosmological hydrodynamical perturbation theory, 
one finds that in all cases of interest, the Bardeen po- 
tential, which represents the typical scalar metric fluctu- 
ations, possesses divergent modes near the bounce that 
cannot be gauged away. For the quantum modes, these 
divergences appear because the so-called "pump" term 
z" I z acts as an infinitely deep potential well : the cor- 
responding D— wave then turns out to be unbounded for 
ry = 0, so that an infinite amount of particles are pro- 
duced. This is similar to the divergent mode in ordinary 
inflationary scenarios, with one fundamental difference : 
in the usual scenarios, the divergence comes from the 
singular point (a = 0) of the metric, and can therefore 
safely be neglected, or altogether ignored, on the basis 
that physics cannot be described by any known theory 
at this point. Near a bounce, however, such an argument 
fails as one may construct perfectly consistent theories 
according to which Einstein gravity and hydrodynamics 
should hold at the bounce. 

The divergent behavior which is obtained in all rele- 
vant cases actually means that any perturbation present 
at the time a bounce occurs would rapidly become non- 
linear. In some cases, this unbounded growth happens 
in the bounce itself. If, in addition, the bounce implies 
a violation of the NEC, non linearities appear in at least 
one other situation : when the NEC is recovered. Hence, 
in many models, one have divergences in two or more mo- 
ments of their history. As such a phenomenon happens 
for all wavelengths (the Hubble radius begin infinite near 
the bounce), one would then expect an ensemble of very 
strong inhomogeneities to be formed prior to any regular 



phase of the Universe. tj 

Note that this does not preclude the actual Universe 
to have a remote bounce. It is possible that after the in- 
homogeneities are formed, a small region of the Universe 
begins to inflate and eventually yields our Universe, much 
like in chaotic inflation, rendering the bounce spacelike 
separated from us (thereby departing from our hypothe- 
sis of physically observable bounce). The bounce in itself 
is very likely to be quite far from any observational veri- 
fication. Its only interest, at the present time, is to reveal 
that it is theoretically possible to construct non singular 
cosmological models. 

Summarizing, a Universe filled only with hydrodynam- 
ical fiuids and with an observable bounce would already 
have been very highly inhomogeneous by the time nucle- 
osynthesis initiates. This conclusion is, from the point 
of view of observational cosmology, inadmissible. This 
shows that one at least of the following statement is true : 

i) No observable bounce ever took place in the Universe, 
or, stated differently, the Universe always expanded. 

ii) If an observable bounce took place, then the dominant 
contribution in the stress-energy sourcing Einstein equa- 
tions was not expressible in the form of a perfect fiuid 
so that hydrodynamic cosmological perturbations do not 
make sense. 

iii) Einstein equations are not valid in the time region 
where the Universe bounced. 

Hence, in the framework of our hypothesis (4- 
dimensional GR and hydrodynamical perturbations), we 
are therefore led to postulate a No-Bounce Conjecture ac- 
cording to which the opening question can be positively 
answered : from a pragmatic point of view, the Universe 
has always expanded. 
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**There is the very specific case where the initial conditions, 
for all wavelength k, are chosen such that the diverging mode 
is vanishing at all times, and if all the solutions of the inhomo- 
geneous equation, when SS ^ 0, are regular. Such an extreme 
choice cannot, of course, be ruled out on purely mathemati- 
cal grounds, but it seems fairly unlikely. In the models where 
the divergences appear more then once in their history, even 
mathematically this may not be possible. 
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APPENDIX : CONNECTING THE BOUNCE 
WITH THE RADIATION DOMINATED ERA OF 
THE STANDARD MODEL 

In this appendix we will prove that (3 < and (3' 
near the bounce when the conditions of case 4 are valid 
{a{ri) = ao + br]'^ + ... in the neighborhood of the bounce). 



the above expansion.tll Physicists believe that no rea- 
sonable dimensionless function f{r]) should contain such 
fine tuned very large a-dimensional numbers^o we con- 
clude that 7(77) = 77 -I- ... for 77 < rjeq ^ lO"!^, and that 
a — \/2bao ~ lO^^cm which implies that 26/ao > 10^^. 
Hence, /3 = 1 — — is negative, and also in this case the 
NEC is violated."" 



To prove that (3 < 0, we will restrict ourselves to the 
single non trivial case IC — 1, where P = 1 — 2b/ a^. 

Considering conformal time ranging from a neighbor- 
hood of the last bounce until standard model evolution (a 
has the minimum value oo in this range), one can define 
without ambiguity the function 



or the inverse relation 



2bao 



a(77) = ^Jal + 2baoP{v) 



(73) 



(74) 



where all initial conditions and physical parameters of 
the solution are contained in b and ao, which means that 
f{ri) is a pure a-dimensional mathematical function. 

Near the bounce we obtain that /(?7) = r] + 0{Tf) 
(remember that this is the case d = Q). Furthermore, 
any such model must be connected with the cosmological 
standard model before nucleosynthesis, which means that 
from rj — r]n ^ 10"'^"'^ to 77 = 77eg ~ 10~^ (equal matter 
and radiation equilibrium) one should have 



0(77) ~ asin(?7) arj 



(75) 



where a ~ 10 cm. Substituting this expression in Eq. 
([73I), and knowing that < a„ ~ 10^^ cm, we get 



fiv) 



\/2ba^ 



(76) 



However, we know that around 77 = the function / is 
exactly 



fiv) + 



m„3 

6! ' 



(77) 



On < a^j and al\/2ha^ > 10^^. 



If 2b /ao < 1 then 2bao < uq 

Knowing that Eq. ( [76[ ) is valid already at 77 ~ 10^^^, the 
compatibility of Eqs. (|7^), ( |7^ ) and 2b /ao < 1 is possi- 
ble only if /'"(O) and higher derivatives be a-dimensional 
huge number, which must be extremely fine tuned in 
order to yield the factor 1 multiplying the t] term in 



We will now prove that (3' ^ around the point where 
(3 changes sign. The change in sign of f3 must occur before 
nucleosynthesis. Hence we can write 



a(?7) = \J2ba0rf 



(78) 



Calculating (3 given in Eq. ( p^ and equating it to zero 
we get the following quartic equation : 

Oa^b^^^ ^ ^i)al{K,ao + 26)77^ -|- al{]Cao - 26) = . (79) 

If /3' = when /3 = 0, then /3" must also be zero. If it is 
not, then (3 is not changing sign. Hence, (3 must be of the 
form (3 = A{ri — ?7i)(77 — 770)'^ if (3' ~ 0, where 770 is some 
root of the quartic equation (fz^), r/i ^ 770, and ^ is a 
constant. But the polynomial form of Eq. (^9|) does not 
contain any cubic 77^^ term so that (3 cannot have three 
equal roots. Hence, we must conclude that ^ 0. 
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